Electronic structure, exchange interactions and Curie temperature in diluted III-V 
magnetic semiconductors: (GaCr)As, (GaMn)As, (GaFe)As. 



L.M. Sandratskii and P. Bruno 
Max-Planck Institut fur Mikrostrukturphysik, D-06120 Halle, Germany 

We complete our earlier (Phys. Rev. B, 66, 134435 (2002)) study of the electronic structure, 
exchange interactions and Curie temperature in (GaMn)As and extend the study to two other di- 
luted magnetic semiconductors (GaCr)As and (GaFe)As. Four concentrations of the 3d impurities 
are studied: 25%, 12.5%, 6.25%, 3.125%. (GaCr)As and (GaMn)As are found to possess a num- 
ber of similar features. Both are semi-metallic and ferromagnetic, with similar properties of the 
interatomic exchange interactions and the same scale of the Curie temperature. In both systems 
the presence of the charge carriers is crucial for establishing the ferromagnetic order. An important 
difference between two systems is in the character of the dependence on the variation of the number 
of carriers. The ferromagnetism in (GaMn)As is found to be very sensitive to the presence of the 
donor defects, like AsGa antisites. On the other hand, the Curie temperature of (GaCr)As depends 
rather weakly on the presence of this type of defects but decreases strongly with decreasing number 
of electrons. We find the exchange interactions between 3d atoms that make a major contribution 
into the ferromagnetism of (GaCr)As and (GaMn)As and propose an exchange path responsible 
for these interactions. The properties of (GaFe)As are found to differ crucially from the proper- 
ties of (GaCr)As and (GaMn)As. (GaFe)As does not show a trend to ferromagnetism and is not 
half-metallic that makes this system unsuitable for the use in spintronic semiconductor devices. 



I. INTRODUCTION 

The discovery of the ferromagnetism in (GaMn)As 1 
with the Curie temperature of Gao.947Mno.o53As as high 
as 110 K attracted much attention to the III-V diluted 
magnetic semiconductors. This attention was stimu- 
lated by both the possibility of deeper understanding of 
the fundamentals of long-range ferromagnetism in semi- 
conductors and the practical need for the ferromagnetic 
diluted magnetic semiconductor (DMS) materials with 
high Curie temperature for the realization of the semi- 
conductor spin-electronic devices. 

On the theoretical side both the model-Hamiltonian 
approach (see recent reviews^) and the parameter-free 
calculations within the density functional theory (see, 
e.g., a review^) are intensively used to study the mag- 
netism of the DMS. One of the important directions of 
the both theoretical techniques is the calculation of the 
exchange interactions and Curie temperature in the DMS 
systems (see, e.g.^SiS*2iili). 

In the previous paper 10 (hereafter referred to as I) 
we used the supercell and frozen-magnon approaches to 
study the exchange interactions and Curie temperature 
of (GaMn)As for four different Mn concentrations. The 
present work extends the study to the cases of Cr and Fe 
impurities. Since a Cr atom has one electron less and an 
Fe atom one electron more than a Mn atom, comparison 
of the results for three DMS systems allows the investiga- 
tion of the trends in the variation of magnetic properties 
within the series of the DMS. 

One of the issues attracting much attention is the de- 
pendence of the Curie temperature on the presence of 
nonmagnetic defects. Such defects can be introduced 
purposely, e.g., by co-doping, or appear in an uncon- 
trolled manner during the sample preparation (e.g., AsGa 
antisites). Since these defects change the density of car- 



riers they can strongly influence the magnetism of the 
DMS systemsiSii4 This influence is also addressed in the 
paper. 

The Cr and Fe impurities in GaAs were recently stud- 
ied theoreticallyiSiii In RefJi the exchange interactions 
within small clusters of the 3d impurities have been 
studied. In Ref^i, the coherent potential approxima- 
tion has been used to estimate the energy difference be- 
tween the ferromagnetic and spin-glass states. In both 
works the authors came to the conclusion of the strong 
ferromagnetic exchange interactions in (GaCr)As. No 
trend to the ferromagnetism has been found in (GaFe)As. 
First attempts to synthesize these systems have been 
rcported>i2i±£ Although the dominant magnetic interac- 
tions between Cr atoms were found to be ferromagnetic 
no high Curie temperature in (GaCr)As was detected^ 
Interestingly, the zinc-blende CrAs that has recently been 
successfully grown on GaAs is ferromagnetic with the 
Curie temperature higher than room temperature^ This 
system has been designed on the basis of the calcula- 
tions within the framework of the density-functional the- 
ory (DFT) and its properties are is good agreement with 
theoretical predictions^ 

The purpose of the present work is to provide detailed 
study of the exchange interactions and Curie tempera- 
ture in the series of diluted magnetic semiconductors: 
(GaCr)As, (GaMn)As, (GaFe)As. Much attention is de- 
voted to the comparative analysis of the systems. One of 
our aims is to further stimulate experimental interest to 
(GaCr)As and (GaFe)As systems. 

The remainder of the paper is organized as follows. In 
Sect. II we discuss a simple two-band model studying 
the relation between the magnetic structure, the band 
occupation and the band energy. In Sect. Ill we briefly 
present the calculational approach. In Sect. IV we dis- 
cuss the calculational results for (GaMn)As, (GaCr)As, 
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and (GaFe)As systems. Our conclusions are given in 
Sect. V. 



II. KINETIC EXCHANGE AND BAND 
OCCUPATION 

It is common to treat the magnetism of the DMS 
in terms of the competition between the antiferromag- 
netic superexchange and ferromagnetic kinetic exchange 
through charge carriers (see, e.g., Refs&i^) The defi- 
nitions of different types of exchange interactions rely 
on different model Hamiltonians and their perturbative 
treatments; 16 Since the DFT is not based on a model 
Hamiltonian approach and does not use a perturbative 
treatment, various exchange interactions appear in the 
calculational results in a mixed form. In this situation 
the studies of simple models of electron systems relevant 
to the problem provide information useful in qualitative 
interpretation of the DFT results. 

Here we consider a simple two-band tight-binding 
model of itinerant electrons experiencing local exchange 
fields of atomic magnetic moments. We consider helical 



configurations of the atomic moments and study the de- 
pendence of the band energy of the system on the mag- 
netic structure. Three cases are discussed: completely 
filled bands, almost empty bands and almost filled bands. 
The helical structures are defined by the formula 

e„ = (cos(q • R„) sin#, sin(q • R„) sin 6*, cos6>) (1) 

where R„ are the lattice vectors, q is the wave vector of 
the helix, e ra is the unit vectors in the direction of the 
magnetic moment at site R„, polar angle 9 gives the de- 
viation of the moments from the z axis. The helical struc- 
tures allow to describe broad range of magnetic configu- 
rations from collinear ferromagnetism (9 — or q = 0) 
to collinear antiferromagnetism (q = and 9 = 90°, 
K is a reciprocal lattice vector). 

The tight binding method for spiral structures was dis- 
cussed in its general form in Reftf . By neglecting the 
difference in the spatial dependence of the basis func- 
tions with opposite spin projections and by preserving 
only the single-center matrix elements of the exchange 
potential we arrive at the following simple form of the 
secular matrix: 



cos 2 (f )H- + sin 2 (f )#+ - f -± sin6>(iJ_ - H+) 



■±sin0(i2_ - ff 4 



wafmH- +cos 2 (f)ii 4 



A 
2 



(2) 



where #_ = H(k- |q), H+ = H(k+ ±q), and H(k) de- 
scribes spin-degenerate bands of a non-magnetic crystal. 
The derivation of Eq. J5J is based on the special symme- 
try of the helical structures^ This model describes one 
electron band with two opposite spin projections. At 
each lattice site R n the band-electron states experience 
a local exchange field of strength A. The model neglects 
the effects of hybridization between magnetic and band 
electrons and therefore is not suitable for the descrip- 
tion of the superexchange^ It, however, takes into ac- 
count interactions characteristic for the kinetic exchange 
through charge carriers ^i^S&SL We will study the de- 
pendence of the energy of the band electrons on the mag- 
netic configuration. 

The eigenvalues of the matrix J2J) have the form 

e ± (k) = ±(H_+H + ) (3) 

±{\{H- - H + f - ±Acos9(H_ - H+) + Ja 2 )^ 

and are illustrated in Fig. ^ To understand the struc- 
ture of the electron bands of the spiral it is important to 
notice that the ferromagnetic configuration of the local 
moments can be considered as a spiral with 9 = and 
arbitrary q. The consideration of the ferromagnet as a 
spiral leads to a rigid shift of the spin-up subband to the 



right by 2 and the shift of the spin-down subband to 
the left by ^. If now the angle 9 becomes nonzero the 
shifted bands of the opposite spin hybridize because of 
non-zero off-diagonal elements of the secular matrix [Eq. 
@]. This leads to the repulsion of the states that are 
close in energy (e.g., the states at q = —0.5 in Fig. 
Additionally, there is an intraband mixing of the states 
separated by vector q in the reciprocal space (see the 
diagonal elements of the secular matrix. 

First, let us consider the case of completely filled 
bands. In this case the total energy equals to 

E h = [ dk[(e_(k)+£ + (k)] =2 f dkff(k) (4) 
Jbz Jbz 

and does not depend on the magnetic configuration. (In 
Eq. Q the integration is carried out over the first Bril- 
louin zone (BZ).) This means that the kinetic exchange 
taken into account by Eq. J2J) does not influence the mag- 
netic configuration in the system with filled bands. The 
superexchange not accounted for by Eq.J2J) plays the 
main role in the case of completely filled bands. 

Now let us consider almost empty bands. In this case 
the minimum of the electron energy corresponds to the 
ferromagnetic structure. To prove this we show that any 
deviation from the ferromagnetic structure increases the 
the minimal energy of the band states (see Fig. ^for an 
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FIG. 1: The energy bands (thick solid line) of the spiral struc- 
ture with q = 0.4, A = 2, and 9 = 45°. H(k) = I - cos(Trk). 
The thin solid lines show the bands of a ferromagnetic con- 
figuration. The broken lines give the ferromagnetic bands 
shifted in the reciprocal space according to the given q (see 
the text). The minimal energy of the ferromagnetic config- 
uration is lower than for the spiral. On the other hand, the 
maximal energy of the ferromagnet is higher than for the spi- 
ral. 




illustration) . 

It is sufficient to consider the lower branch of the eigen- 
values £-(k). Let us take an arbitrary k. If If ff_(k) < 
H+(k) than g_(k) > ff-(k) - ±A. If i?_(k) > H + (k) 
than e_(k) > H + (k) - ±A. Therefore, in both cases the 
eigenstate of the helix £-(k) is higher than one of the 
states of the ferromagnetic crystal. For the k vector sat- 
isfying H_(k) = H + (k) = h, e_(k) = h - |A. In the 
usual case of H (k) having minimum at the center or the 
boundary of the BZ, h — ^A is not the minimal energy 
of the ferromagnetic structure. 

Also for almost filled bands the configuration of low- 
est energy is ferromagnetic since any deviation from the 
ferromagnetic state increases the energy of the electrons. 
To show this we first notice that the maximal energy of 
the band states always decreases with deviation from fer- 
romagnetism. The proof of this is similar to that given 
above for the minimal energy. Since the energy of the 
completely filled bands does not depend on the magnetic 
configuration, the highest energy of the hole states means 
the lowest energy of the system. 

In the case of exchange splitting A much bigger than 
the width of the nonmagnetic band H (k) the low-energy 
eigenvalue of the secular matrix given by Eq. © takes the 



form 

e_(k) = cos 2 (^)[ff(k-iq)-^] (5) 

+ s in 2 (~)[ff(k+iq)-|] 

It can be easily verified that this subband fulfills all three 
properties formulated above for the set of two bands. 

These properties of the model agree qualitatively with 
those obtained within other models, e.g., within the 
model of two interacting impurities^ We will refer to 
these properties in the discussion of the results of the 
DFT calculations in the following sections. 

III. CALCULATIONAL APPROACH 

In the calculations we use the scheme discussed in I. 
This scheme is based on the supersell approach where 
one of the Ga atoms in a supercell of zinc-blende GaAs is 
replaced by the 3d atom. The calculations are performed 
for four values of the concentration x: 0.25, 0.125, 0.0625, 
and 0.03125. 

The calculations were carried out with the augmented 
spherical waves 2 ^ (ASW) method within the local density 
approximation (LDA) to the DFT. In all calculations the 
lattice parameter was chosen to be equal to the experi- 
mental lattice parameter of GaAs. Two empty spheres 
per formula unit have been used in the calculations. The 
positions of empty spheres are (0.5, 0.5, 0.5) and (0.75, 
0.75, 0.75). Radii of all atomic spheres were chosen to 
be equal. Depending on the concentration of Mn, the 
super cell is cubic (x=25%, a x a x a, and x=3.125%, 
2a x 2a x 2a) or tetragonal (x=12.5%, a x a x 2a and 
6.25%, 2a x 2a x a). 

The densities of states (DOS) presented in the paper 
are calculated for self-consistently determined ferromag- 
netic states of the corresponding systems. 

To describe the exchange interactions in the system we 
use an effective Heisenberg Hamiltonian of classical spins 

H eff = -^Jijti-ej (6) 

where Jy is an exchange interaction between two 3d sites 
and ej is the unit vector pointing in the direction 
of the magnetic moments at site i. 

To estimate the parameters of the Mn-Mn exchange 
interaction we performed calculation for the following 
frozen-magnon configurations: 

6i = const, fa = q • R, (7) 

where 9i and 4>i are the polar and azimuthal angles of 
vector e», R^ is the position of the ith Mn atom. The 
directions of the induced moments in the atomic spheres 
of Ga and As and in the empty spheres were kept to be 
parallel to the z axis. 
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It can be shown that within the Heisenberg model © 
the energy of such configurations can be represented in 
the form 



E(d,<i)=Eo(6)~9 2 J(<i) 



(8) 



where Eq does not depend on q and J(q) is the Fourier 
transform of the parameters of the exchange interaction 
between pairs of Mn atoms: 



J(q) = J oj cx P(«q ' R oj)- 



(9) 



In Eq. JSJ angle is assumed to be small. Using J(q) one 
can estimate the energies of the spin- wave excitations: 



w (q) = — [J(0)-J(q)l = , 2 r 

M WJ M(l-cos0) 



(£(0,q)-£J(0,O)) 
(10) 

where M is the atomic moment of the Mn atom. Per- 
forming back Fourier transformation we obtain the pa- 
rameters of the exchange interaction between Mn atoms: 



(11) 



FIG. 2: Tc of (Ga,Mn)As with different Mn concentrations 
as a function of the electron number n. n — corresponds to 
the system Gai-^Mnz As with no additional donor or acceptor 
defects. 
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The calculation oiE(9, q) for different Mn concentrations 
has been performed for uniform meshes in the first BZ. 
Angle is selected in the proportionality region between 
[£(0,q) -£(0,0)] and (1 -cos0). 

The Curie temperature was estimated in the mean-field 
(MF) approximation 



k B T$ p 



(12) 



We use rigid band approach to calculate the exchange 
parameters and Curie temperature for different electron 
occupations in each of the systems studied. This allows 
to simulate the influence of the As antisites and the non- 
magnetic co-doping for a given concentration of 3d im- 
purities. We assume that the electron structure calcu- 
lated for a DMS with a given concentration of the 3d 
impurity is basically preserved in the presence of defects. 
The main difference is in the occupation of the bands 
and, respectively, in the position of the Fermi level. For 
each electron occupation we calculated the energy of the 
frozen-magnon states and the interatomic exchange pa- 
rameters. 



IV. CALCULATIONAL RESULTS 
A. (GaMn)As 

The main body of the calculational results for 
(GaMn)As is presented in I. Here we discuss the depen- 
dence of the calculated Curie temperature on the num- 
ber of carriers. The experimental studies show that the 



concentration of holes in (Ga,Mn)As is lower than the 
concentration of the Mn atomsS One of the important 
factors leading to the low concentration of holes is the 
presence of the As antisites. (See Reism^^^^ for ear- 
lier DFT studies of the influence of the As antisites on 
the properties of the DMS.) 

Since As has two more valence electrons compared with 
Ga, each As antisite compensate the holes produced by 
two Mn atoms. We study the influence of the number of 
electrons in the broad interval from n = —2 (two elec- 
trons per supercell less) to n = 2 (to electrons per super- 
ccll more). Concretely, for each frozen-magnon config- 
uration studied we perform the total energy calculation 
for different n and then for each n the procedure of the 
evaluation of the exchange interactions described in Sec. 
IIIII is used. Negative values of the Curie temperature in 
Fig|5] indicate an instability of the ferromagnetic state 
due to dominating antiferromagnetic interactions. 

The results of the calculations are shown2& in Fig. [21 
The following features are important. First, for all con- 
centrations x the value of the Curie temperature oscil- 
lates with occupation. It is negative for n — —2 (two 
additional holes), increases with increasing occupation, 
changes sign and reaches the maximum in the interval 
between n — — 1 (one additional hole) and n = 0. Then 
it decreases again, changes sign in the interval between 
ri = and n = 1 (one additional electron) reaches mini- 
mum and increases again. The only deviation from this 
scenario is the case of x — 3.125% where the curve has 
two maxima, one close to n = 1 and another to n = 0. 
(The two-maxima shape of the curve can be related to 
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the properties of the electron DOS presented in I in the 
energy region just below the Fermi level: the closeness of 
the top of the spin-down valence band to the Fermi level 
and a deep minimum in the spin-up DOS.) A remark- 
able feature in the dependences is the kink at n = 1 for 
x = 3.125% and x = 6.25%. This n value corresponds 
to a completely filled valence band and empty conduc- 
tion band. Correspondingly, for n < 1 the carriers are 
holes in the valence band and for n > 1 the carriers are 
electrons in the conduction band. The abrupt change of 
the character of the carrier states at n = 1 results in the 
discontinuity of the slope of Tc{n). For x — 12.5% and 
x = 25% the kink is not obtained because of the overlap 
of the valence and conduction bands. Also the minimum 
for these concentrations is shifted to a non-integer value 
of n that depends on the details of the band overlap and 
cannot be predicted without calculations. 

The amplitude of the curves decreases with decreasing 
x. Therefore, potentially the highest Curie temperature 
can be reached for the highest Mn concentration. This, 
however, needs large-scale tuning of the number of car- 
riers. (Both, the growth of the samples with Mn con- 
centration of the order of 20% and corresponding large- 
scale tuning of the number of carriers are technologically 
hardly possible.) It is interesting to consider the depen- 
dence of Tq on concentration x for various numbers of 
electrons n. Since the form of the curves Tc(n) differs 
considerably for different x they intersect (Fig. [5J. The 
intersection of the curves influences the dependence of 
Tc on concentration x (Fig. [31 . Indeed, if before inter- 
section at some n a , Tc for concentration x\ is larger than 
Tc for concentration X2 after the intersection the relation 
is opposite. The last property is reflected in Fig. [3] where 
the value of T C AF as a function of x is shown for different 
hole concentrations. Indeed, a monotonous dependence 
for n=— 0.4 is replaced by a nonmonotonous behavior for 
n = with the maximum at x = 12.5%. For n = 0.4 and 
n = 0.6 the maximum corresponds to x — 6.25%. The 
experimental Tc values obtained in Ref^£ are in good 
agreement with the values for n = 0.6 . For n = 0.8 the 
maximum value of the Curie temperature corresponds to 
the lowest Mn concentration. 

Note that the mean-field approximation usually over- 
estimates the Curie temperature and the comparison of 
Tq F with experiment should be made with caution. 
For example, the value of the Curie temperature of 
(GaMn)As obtained in I within the random phase ap- 
proximation (RPA) is about 15-20% smaller than the 
T(f F . This means that good agreement with the ex- 
perimental Curie temperatures just for n — 0.6 might be 
accidental. The most important result of the calculation 
is the qualitative trends obtained with the variation of 
the number of carriers. 

In general, Figs. 12131 confirm the picture of the ferro- 
magnetism in DMS as mediated by the charge carriers. 
Indeed, the decrease of the number of holes in the va- 
lence band (n > 0) leads to a fast decrease in the trend 
toward ferromagnetism. However, because of the oscil- 



FIG. 3: Tc of (Ga,Mn)As as a function of the Mn concen- 
tration x for different electron numbers n. n is defined as in 
Fig. |5| Note that the number of holes p per Mn atom equals 
1-n. The stars show the experimental values of the Curie 
temperature. 
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lating character of the curves Tc(n), the increase of the 
number of holes to more than one per Mn (n < 0) has 
only small potential for obtaining higher Curie tempera- 
ture in (GaMn)As. 

Figure|2(and the corresponding Figs. H0landlT3lin the 

following) show that the contribution of the completely 
filled bands into exchange interaction is always antiferro- 
magnetic. Indeed, in all cases of completely filled bands 
(n = 1 in Fig. , n = — 1 and n = 2 in Fig. E3 0r n = 2 
in Fig. I13JI the estimation of the Tq F is negative. On 
the other hand, the deviation from the completely filled 
bands to smaller n (creating holes in the filled bands) or 
to larger n (creating electrons in the empty bands) leads 
to increasing Tc IF and, respectively, to increasing trend 
to ferromagnetism. The results of the study of the two- 
band model in Sect. [n]hclp to understand the physical 
mechanism of this property in terms of the changes in the 
band structure with the deviation of the magnetic con- 
figuration from the collinear ferromagnetic one. Indeed, 
we have seen that for completely filled bands the changes 
in the energy of different states compensate and no in- 
crease of the total energy is connected with the deviation 
from the collinear ferromagnetic structure. The antifer- 
romagnetic superexchange obtained in Fig. [21 for n = 1 
is not described by the model of Sect. [H] since the hy- 
bridization between magnetic and band electrons is not 
taken into account. For the case there are partially field 
band with holes or electron present, the ferromagnetic 
state becomes energetically preferable. These properties 
correlate with the view of the magnetism of the DMS 
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FIG. 4: The DOS of Gai-^Cr^As. The DOS is given per 
unit cell of the zinc-blende crystal structure. The DOS 
above(below) the abscissas axis corresponds to the spin- 
uptown) states. 
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as governed by the competition between the antiferro- 
magnetic exchange through completely filled bands and 
ferromagnetic exchange through charge carriers (holes or 
electrons). 

B. (GaCr)As 

In Fig. Hwe show the DOS of (GaCr)As for four val- 
ues of the Cr concentration. There is substantial differ- 
ence between these DOS and the DOS of (GaMn)As pre- 
sented in I. In particular, an important different feature 
is a higher energy position of the Cr 3d states relative 
to the GaAs states if compared with the position of the 
Mn 3d states in (GaMn)As. For the Cr concentration 
of x = 3.125% this energy shift results in an impurity 
band lying within the semiconducting gap of GaAs and 
separated in energy from both the valence and the con- 
duction bands. (In the case of (GaMn)As corresponding 
states lie at the top of the valence band.) The impurity 
band has spin-up character. The replacement of one Ga 
atom in the supercell of GaAs by a Cr atom does not 
change the number of the spin-down states in the va- 
lence band. In the spin-up channel there are, however, 



TABLE I: Magnetic moments in Gai-zCr^As and 
Ga,i- X ¥e x As. There are shown the moment on the 3d 
impurity, the induced moment on the nearest As atoms, and 
the magnetic moment of the supercell. All moments are in 
units of /is- 







0.25 


0.125 


X 

0.0625 


0.03125 




Cr 


3.10 


3.11 


3.13 


3.15 


(GaCr)As 


As 


-0.09 


-0.08 


-0.07 


-0.07 




cell 


3.00 


3.00 


3.00 


3.00 




Fe 


3.18 


3.29 


3.40 


3.46 


(GaFe)As 


As 


0.03 


0.04 


0.05 


0.06 




cell 


3.61 


3.87 


4.15 


4.39 



five additional energy bands which are related to the Cr 
3d states. Two of them lie within the valence band and 
three in the semiconducting energy gap. Since there are 
five extra energy bands and only three extra electrons 
(the atomic configurations of Ga and Cr are 4s 2 4p 1 and 
3d 4 4s 2 ) the impurity band is not filled. The integrated 
number of the occupied (electron) states in the impurity 
band is one per Cr atom (correspondingly, one electron 
per supercell). On the other hand, the integrated num- 
ber of empty (hole) states in the impurity band is two 
per Cr atom (two holes per supercell). 

With increasing x, the impurity band becomes 
broader. At x = 6.25% it touches the spin-down states of 
the conduction band still being separated by an energy 
gap from the valence band. At x = 12.5% and x = 25% 
the impurity band overlaps with both valence and con- 
duction bands. For all concentrations the Fermi level lies 
in the energy gap of the spin-down DOS. Therefore for 
all x the system is half-metallic and is characterized by 
a 100% spin-polarization of the carriers at the Fermi en- 
ergy. The latter property is very important for efficient 
spin-injection to semiconductors. 

In Table [I] we present the magnetic moment of the Cr 
atom, the induced moment on the nearest As atom, and 
the magnetic moment per supercell. Both the Cr and 
As moments vary weakly with the change of x. This 
is the result of the half-metalicity that leads to a fixed 
number of the spin-up and spin-down electrons in the 
system and, correspondingly, complicates the change of 
atomic moments. 

Taking as an example the system with x — 3.125% we 
find that the contribution of the 3d electrons into the Cr 
moment is 3.10/is with the rest 0.05/is coming from the 
4s and 4p electrons. The total number of the 3d electrons 
in the Cr sphere is 4.38. The main contribution to the 
induced As moment comes from the 4p electrons. 

Interestingly, the induced As moment in (GaCr)As 
is larger than in (GaMn)As though the Mn moment is 
larger than the Cr moment. The explanation of this fact 
is in the hybridization of the spin-up states of the 3d 
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FIG. 5: Distribution of the carrier states over atoms (up- 
per panel) and induced magnetic moments (lower panel) 
for Gai-^Crj; As and Gai_ I Mn I As with x = 3.125%. For 
(GaCr)As, the distribution is given for the occupied (elec- 
tron) states in the impurity band. For (GaMn)As,the distri- 
bution is given for empty (hole) states in the valence band. 
All values are for the coordination spheres. The numbers of 
atoms in the coordination spheres are shown at the top of the 
picture. The calculated values for (GaMn)As are taken from 
I. 
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and As atoms. Because of this hybridization, part of the 
As spin-up states contributes to the hole states in the 
valence (GaMnAs) or impurity (GaCrAs) band and is 
unoccupied. The number of hole states in (GaCr)As is 
larger than in (GaMn)As, therefore the negative induced 
moment on the As atom is larger as well. 

Since the ferromagnetism is expected to be mediated 
by carriers it is important to consider their localization 
about Cr atoms. In Fig. [5] we show for x = 3.125% the 
distribution in the supercell of the one electron in the oc- 
cupied part of the impurity band. This distribution was 
calculated as follows. At first the distribution over atoms 
was found for each state of the occupied part of the im- 
purity band. Then the summation over these states was 
performed. The number of electrons in the occupied part 
of the impurity band is exactly one per Cr atom (or, what 
is equivalent, per supercell). For comparison, we plot in 



Fig. El the data for one hole per Mn atom in (GaMn)As 
calculated similarly. About 44% of the electron is in the 
Cr sphere. The first coordination sphere of As contains 
13% of the state. The third coordination As sphere with 
distance 1.09a from the Cr atom contains 11% of the elec- 
tron. We will show that the large portion of the electron 
states at the atoms of the third coordination sphere of 
the As atom is important for the magnetic ordering in 
the system. The first coordination sphere of Ga atoms 
contains 9% of the state. Each of the other coordination 
spheres contains less that 3% of the state. 

Another important characteristic for establishing the 
exchange interaction between the Cr atoms is the value 
of induced moments (FiglSJ. The induced moments on all 
As atoms are antiparallel to the Cr moment. The same 
feature was obtained in I for (GaMn)As (see Fig|S} and is 
explained by the contribution of the As spin-up states to 
the hole states in the impurity band discussed above. On 
the Ga atoms, the induced moments are parallel to the Cr 
moment. Note that the parts of the states corresponding 
to the Ga atoms and the induced moments on the Ga 
sites are very similar in both (GaCr)As and (GaMn)As. 
On the other hand, for the 3d and As atoms there is 
substantial difference between two systems. 

For (GaCr)As, there is a clear correlation between the 
value of the induced moment and the part of the elec- 
tron states in the impurity band corresponding to a given 
atom (FiglSJ). In particular there is a sizable induced mo- 
ment at the third coordination sphere of As atoms. Note 
that the hole states in (GaMn)As are less localized about 
Mn atom compared with the electron states localization 
about Cr atoms in (GaCr)As. On the other hand, the 
values of the induced moments are larger in (GaCr)As. 
Thus apriori it is not clear which of the two systems has 
potential for the ferromagnetism with higher Curie tem- 
perature. 

In Fig. we show calculated interatomic exchange 
parameters. The quantitative comparison with the ex- 
change parameters of (GaMn)As presented in I shows 
that they are substantially different. There is however 
a number of qualitative similarities. First, in both cases 
the Heisenberg model with the interaction between the 
first nearest neighbors only is not able to describe the 
magnetism of the system. Second, the exchange interac- 
tions are rather quickly decreasing with increasing dis- 
tance between atoms. In Fig. {7\we plot the dependence 
of the absolute value of the exchange parameters for three 
different impurities on the interatomic distance using log- 
arithmic ordinate axis. The decrease of exchange param- 
eters is close to exponential that is expected for semi- 
conducting and half- metallic systems In Fig. H3 we 
show the variation of the Curie temperature [Eq. l|12fl ] 
with increasing number of the contributing coordination 
spheres. For instance, no noticeable contribution to the 
Curie temperature is obtained from the interactions be- 
tween Mn atoms at the distances larger than 3a. Third, 
the dependence of the exchange parameters on the dis- 
tance between Cr atoms is not monotonous. Forth, the 
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FIG. 6: The parameters of the exchange interaction be- 
tween Cr atoms (upper panel) and the variation of Tq F with 
increasing number of the contributing coordination spheres 
(lower panel). The abscissa gives the radius of the coordina- 
tion sphere in the units of the lattice parameter of the zinc- 
blende crystal structure. The calculations has been performed 
up to r = 7a. Because of very small values the exchange 
parameters for large r are not shown. The broken vertical 
lines show the coordination spheres with strongest exchange 
interaction with the central atom. Two data points for the 
same system and the same radius (e.g., at r = 2.828a for 
x = 12.5% and ) mean that there are two inequivalent coor- 
dination spheres with equal radii. 
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coordination spheres which contribute importantly into 
the magnetism of (GaMn)As provide, in most cases, such 
contribution also in (GaCr)As. 

In Fig. [S] the Curie temperature of (GaCr)As for 
four values of the Cr concentration is presented. For 
x = 3.125% and x = 25% the Curie temperature is higher 
than in (GaMn)As, for x = 6.25% and x = 12.5% the 
relation is opposite. In general both systems give Tq 
in the same temperature interval. Only for very large 
x the Curie temperature of (GaCr)As becomes substan- 
tially larger. 

The calculated exchange parameters combined with 
the information about the induced moments and spatial 
distribution of the states of the impurity band allow to 
suggest the exchange path responsible for the realization 
of the ferromagnetic state. 



FIG. 7: The absolute values of the exchange parameters 
between 3d impurities as a function of the interatomic dis- 
tance. The impurity concentration is 25%. The [110] crystal- 
lographic direction is presented. A logarithmic ordinate axis 
is used to visualize fast decay of the exchange parameters. 
The small value of the exchange parameter of (GaFe)As at 
r = 2\/2a results from the oscillation of the parameter as a 
function of r(Fig. I12H . The dashed straight line is a guide for 
the eye. 
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The analysis of the interatomic exchange interactions 
and their contribution into the Curie temperature (Fig. 
|BJ) shows that the main role for all concentrations of 
the Cr atoms is played by the exchange interaction be- 
tween the Cr atoms separated by the vectors parallel to 
(1,1,0) (or to the vectors crystallographically equivalent 
to (1,1,0)). Indeed, for concentrations x = 25%, 12.5% 
the main contribution to the Curie temperature comes 
from the atoms separated by (1,1,0). For lower concen- 
trations the main role is played by the interaction be- 
tween atoms separated by (2,2,0). 

This result together with the information about the in- 
duced moments and the localization of the carrier states 
allow to propose the exchange path for the realization of 
the strong exchange interactions. Let us consider again 
the case of x = 3.125% (Fig. |5J). On the straight line 
connecting the Cr atoms at (0,0,0) and (2,2,0) there are 
three Ga atoms at (±, ±,0), (1, 1,0), and (§, |,0). (The 
first and the third atoms are equivalent.) Since, accord- 
ing to Fig. El only very small portion of the carrier states 
corresponds to these atoms and the induced moments on 
these atoms are also very small they cannot efficiently 
mediate the exchange interaction between Cr atoms. 

On the other hand a substantial part of the elec- 
tron states in the impurity band is on the As atoms 
at and (|, |, |). Simultaneously, these atoms 

carry sizable induced moments. On the basis of 
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FIG. 8: Tq IF of (Ga,Cr)As, (Ga,Mn)As, and (Ga,Fe)As for 
different concentrations of the 3d impurity. 



FIG. 10: Tc of (Ga,Cr)As with different Cr concentrations 
as a function of the electron number n. n = corresponds to 
the system Gai-ajCr^As with no additional donor or acceptor 
defects. 
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FIG. 9: The z — and z = j planes of the supercell for the 
concentration of magnetic impurity of 3.125%. The exchange 
path shown by the thick line connecting the impurity atoms 
goes over As atoms. 
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this information we suggest the following exchange 
path: Cr(0,0,0)-As(i, ± \) - As(f, f , ±)-As(f , f , \) 
-As(|,|,i)-Cr(2,2,0). Here As(|,|,|) is equivalent 
to As(±, \, \) and As(§, f , i) is equivalent to As(|, |, i). 
The calculational information for (GaMn)As (see I) al- 
lows to suggest the same exchange path also for this sys- 
tem. 

To get deeper insight into the magnetism of (GaCr)As 



we calculated the dependence of the Curie temperature 
on the electron occupation (Fig. I10[l . Many qualita- 
tive features of these dependences are similar to those 
for (GaMn)As (Fig. The dependences show oscil- 

lating behavior with the amplitude decreasing with de- 
creasing Cr concentration. For x — 3.125,6.25,12.5% 
there is a kink at n = — 1 similar to the kink at n = 1 
for (GaMn)As. The origin of the kink is also similar to 
that for (GaMn)As and consists in the transition from 
the occupation of the valence-band states at n < — 1 to 
the occupation of the impurity-band states at n > — 1. 
At n — — 1 the valence band is completely filled and the 
impurity band is empty. The curves for x — 3.125% and 
x = 6.25% have another kink at n = 2 which marks the 
end of the filling of the impurity band and beginning of 
the filling of the conduction band. The reason for the in- 
crease of the interatomic ferromagnetic interactions with 
appearance of holes in valence band or electrons in con- 
duction band was discussed in Sect. [H]in terms of two- 
band model and consists in the increase of the energy 
of the system with deviation of the moments from the 
parallel directions. 

An important conclusion from the comparison of the 
Figs. El and [21 is that the role played in the case of 
(GaMn)As by holes in the valence band is played in 
(GaCr)As by the electrons in the impurity band. In par- 
ticular, the decrease of the number of electrons results 
in fast decrease of the Curie temperature similar to the 
behavior obtained in (GaMn)As with decreasing number 
of holes. On the other hand, the decrease of the number 
of holes in (GaCr)As (the region of positive n in Fig ll(J|) 
first leads to a moderate increase of the Curie tempera- 
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FIG. 11: The DOS of Gai-^Fe^As. The DOS is given 
per unit cell of the zinc-blende crystal structure. The DOS 
above(below) the abscissas axis corresponds to the spin- 
uptown) states. 
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ture. Only for larger n (n > 1 for x = 3.125, 6.25, 12.5%) 
a strong decrease of Tc is observed. (Again, a similar 
behavior is obtained for (GaMn)As but with n decreas- 
ing lower than —1.) Therefore the Curie temperature in 
(GaCr)As is less sensitive to the presence of the As an- 
tisites and other nonmagnetic donor defects. Since the 
formation of the As antisites is difficult to avoid this prop- 
erty is of technological importance. 



C. (GaFe)As 

An Fe atom has one additional valence electron com- 
pared with a Mn atom. Thus no hole or electron carriers 
can be expected in (GaFe)As on the basis of the elec- 
tron count. The system can be expected to be strongly 
antifcrromagnetic. However, the calculations show that 
the electron structure of (GaFe)As differs strongly from 
the electron structure of both (GaMn)As and (GaCr)As 
and a simple consideration on the basis of the electron 
count does not apply. The origin of this difference is in 
the lower energy position of the Fe 3d states with respect 
to the GaAs states if compared with the position of the 
Mn states (Fig. 111(1 . As a result, the spin-down impu- 
rity states lie at the top of the valence band and become 
partly occupied. Correspondingly, the spin-up impurity 
states at the top of the valence band are not completely 
filled. The presence of both spin- up and spin-down states 
at the Fermi level is disadvantageous for the realization 
of an efficient spin-injection that needs very high spin- 
polarization of the states at the Fermi level. With in- 



creasing x the impurity bands broaden. For none of the 
concentrations the state of the system is half-metallic. 

In Table H| the Fe magnetic moment as well as the in- 
duced moment on the first As atom are shown. Com- 
pared with (GaCr)As, the moments strongly vary with 
concentration. This variation is a consequence of the 
presence of both spin-up and spin-down states at the 
Fermi level that leads to strong electron transfer between 
two spin channels. The induced moment is now positive 
and its value increases with decreasing x. The reason for 
the change of the sign of the induced As moment is di- 
rectly related to the increased number of 3d electrons in 
an Fe atom compared to Mn and Cr atoms. Indeed, in 
(GaFe)As the number of the hole spin-up states becomes 
very small. These hole states play a crucial role in the 
formation of the negative induced moments in (GaMn) As 
and (GaCr)As since they decrease the number of the As 
spin-up electron. Simultaneously, the number of the As 
spin-down electrons in (GaFe)As decreases since, because 
of a lower energy position of the Fe 3d spin-down states, 
the As spin-down states hybridize with empty impurity 
states. 

The calculation of the exchange interactions shows 
(Fig. I12fl that they are much more antifcrromagnetic 
than for (GaMn)As and (GaCr)As. Correspondingly, the 
estimation of the Curie energy (Fig. [SJ is negative for 
x = 3.125% and x = 25% and close to zero for x = 12.5%. 
It is positive for x = 6.25% but much smaller than the 
Curie temperature for the same concentration of the 3d 
element in (GaMn)As and (GaFe)As. The explanation 
for such behavior is in small number of the spin- up holes. 
The main number of the states at and above the Fermi 
level are of the spin-down type. The spin down hole states 
are strongly localized on the Fe atom (76.5%) and do 
not mediate ferromagnetic interactions efficiently. The 
curves of Tc{n) (Fig. I13fl show that there is a general 
trend to further strengthening of the antifcrromagnetic 
interactions with increasing n. In this case the valence 
band is filled at n — 2. In agreement with the two-band 
model of Sect.[n]the curves have minimum at this point. 
Appearance of holes (n < 2) leads to increasing ferromag- 
netic interactions. Although strong acceptor doping has 
some potential for increasing ferromagnetic interactions, 
(GaFe)As is not a good candidate for the application in 
the semiconductor spin-electronic devices. 



V. CONCLUSIONS 

In the given paper we complete our study of the elec- 
tronic structure, exchange interactions and Curie tem- 
perature of (GaMn)As started in I and report calcula- 
tions for (GaCr)As and (GaFe)As. A different number of 
the valence electrons in Cr, Mn and Fe allows the inves- 
tigation of the trends within the series. We find that 
there is a similarity in the character of the exchange 
interactions of (GaMn)As and (GaCr)As. The depen- 
dence of the exchange interactions and Curie temper- 
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FIG. 12: The parameters of the exchange interaction be- 
tween Fe atoms (upper panel) and the variation of Tq IF with 
increasing number of the contributing coordination spheres 
(lower panel). The abscissa gives the radius of the coordina- 
tion sphere in the units of the lattice parameter of the zinc- 
blende crystal structure. The decay of the exchange parame- 
ters is illustrated in the logarithmic scale in Fig. |7] 
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temperature of (GaCr)As depends rather weakly on the 
presence of this type of defects but decreases strongly 
with decreasing number of electrons. The properties of 
(GaFe)As are found to differ crucially from the properties 

FIG. 13: T^ IF of (Ga,Fe)As with different Fe concentrations 
as a function of the electron number n. n = corresponds to 
the system Gai-^FezAs with no additional donor or acceptor 
defects. 
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ature on the band occupation shows that the presence 
of the charge carriers is crucial for establishing of the 
ferromagnetism in (GaMn)As and (GaCr)As. We find 
that the role played in (GaMn)As by the holes in the 
valence band is played in (GaCr)As by the electrons in 
the impurity band. An important difference between two 
systems is in the character of the dependence on the vari- 
ation of the number of carriers. The ferromagnetism in 
(GaMn)As is very sensitive to the presence of the donor 
defects, like Asc a antisites. On the other hand, the Curie 



of (GaCr)As and (GaMn)As. (GaFe)As does not show 
a trend to ferromagnetism and is not half-metallic that 
makes this system unsuitable for the use in spintronic 
semiconductor devices. 

We show that the sign of the induced moment on the 
As atoms and, correspondingly the sign of the p-d ex- 
change (often denoted as N/3) varies from negative in 
(GaMn)As and (GaCr)As to positive in (GaFe)As. There 
is a connection between the half-metalicity of the system 
and the negative sign of N/3. 

We show that the strongest exchange interactions in 
(GaCr)As and (GaMn)As are between the 3d atoms sep- 
arated by the vector parallel to (1,1,0) or to a vector 
crystallographically equivalent to (1,1,0). On the basis 
of the calculated results we propose the exchange path 
over As atoms that is responsible for the mediating the 
ferromagnetic interactions between the Cr atoms. 
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